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Abstract 

New dark energy models in anisotropic Bianchi type-I (B-I) space-time 
with variable EoS parameter and constant deceleration parameter have 
been investigated in the present paper. The Einstein's field equations 
have been solved by applying a variation law for generalized Hubble's 
parameter in B-I space-time. The variation law for Hubble's parameter 
generates two types of solutions for the average scale factor, one is of 
power-law type and other is of the exponential form. Using these two 
forms, Einstein's field equations are solved separately that correspond to 
expanding singular and non-singular models of the universe respectively. 
The equation of state (EoS) parameter uj is found to be time dependent 
and its existing range for this model is in good agreement with the recent 
observations of SNe la data, SNe la data (with CMBR anisotropy) and 
galaxy clustering statistics. The cosmological constant A is found to be 
a decreasing function of time and it approaches a small positive value at 
late time (i.e. the present epoch) which is corroborated by results from 
recent supernovae la observations. 
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1 Introduction 

Before few years ago two groups (the Supernova Cosmology Project and the 
High-Z Supernova Team) presented evidence that the expansion of the universe 
is accelerating (Garnavich et al. [2H]; Perlmutter et al. [3J 21 [5] ; Riess et al. [5J; 
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Schmidt et al. [7J. These teams have measured the distances to cosmological 
supernovae by using the fact that the intrinsic luminosity of Type la supernovae 
is closely correlated to their decline rate from maximum brightness, which can 
be independently measured. These measurements, combined with red-shift data 
for the supernovae, led to the prediction of an accelerating universe. Both team 
obtained Hm ~ 0.3, Ha ~ 0.7, and strongly ruled out the traditional (Om, ^a) 
— (1,0) universe. This value of the density parameter Ha corresponds to a 
cosmological constant that is small, nevertheless, nonzero and positive, that is, 
A ps 10~ 52 m~ 2 ps 10~ 35 s~ 2 . An intense search is going on, in both theory and 
observations, to unveil the true nature of this acceleration. It is commonly be- 
lieved by the cosmological community that a kind of repulsive force which acts 
as anti-gravity is responsible for gearing up the Universe some 7 billion years 
ago. This hitherto unknown exotic physical entity is termed as dark energy. 

It is an irony of nature and is a puzzling phenomenon that most abun- 
dant form of matter-energy in the universe is most mysterious. The simplest 
Dark Energy (DE) candidate is the cosmological constant A, but it needs to 
be extremely fine-tuned to satisfy the current value of the DE. Alternatively, 
to explain decay of the density, many dynamic models have been suggested, 
where A varies slowly with cosmic time (t) (Overduin and Cooperstock [5]; 
Sahni and Starobinsky [S]; Komatsu et al. [ID])- In addition to models with 
dynamic A(t), many hydrodynamic models with or without dissipative pressure 
have been proposed in which barotropic fluid is the source of DE (Sahni and 
Starobinsky \9 s ; Komatsu et al. [10] ). Chaplygin gas as well as generalized 
Chaplygin gas have also been considered as possible dark energy sources due 
to negative pressure (Srivastava [IT] : Jackiw [12]; Bertolami et al. [13]; Bcnto, 
Bertolami and Sen [T3] ; Bilic, Tupper and Viollier [TS]; Avelino et al. [IS]). 
Other than these approaches, some authors have considered modified gravita- 
tional action by adding a function f(R) (R being the Ricci scalar curvature) 
to Einstein- Hilbert Lagrangian, where f(R) provides a gravitational alternative 
for DE causing late-time acceleration of the universe (Capozziello[I7]; Caroll et 
al. [H] ; Dolgov and Kawasaki [TO] ; Nojiri and Odintsov [30] HJ ; Abdalaa et al. 
[2"2"] ; Mena et al. [53] ). A review on modified gravity as an alternative to DE is 
available in Nojiri and Odintsov|24j. A more comprehensive review is provided 
in Copeland et al. [25]). Recently Gupta and Pradhan [26] have presented an 
entirely new approach that cosmological nuclear energy is a possible candidate 
for DE. All these models are phenomenological in the sense that an idea of DE 
is introduced a priori either in term of gravitational field or non-gravitational 
field. In spite of these attempts, still cosmic acceleration is a challenge. 

At present, there is great interest in cosmological models with a variable 
matter equation of state in the class of equation p — to(t)p (p is the fluid 
pressure and p its energy density). By now, methods allowing for restoration 
of the quantity io(t) from expressional data have been developed (Sahni and 
Starobinsky [27], and an analysis of the experimental data has been conducted 
to determine this parameter as a function of cosmological time (see Sahni et 
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al. [2H] and references therein). DE has been conventionally characterized by 
the equation of state (EoS) parameter mentioned above which is not necessarily 
constant. Recently, the parameter u>(t) is calculated with some reasoning which 
reduced to some simple parametrization of the dependences by some authors 
(Huterer and Turner [25]; Weller and Albrecht [3D]; Chevallier and Polarski [3T] ; 
Krauss et al. 32 : Usmani et al. [33]; Chen et al. [34] . The simplest DE 
candidate is the vacuum energy (to = —1), which is mathematically equivalent 
to the cosmological constant (A). The other conventional alternatives, which 
can be described by minimally coupled scalar fields, are quintessence (w > — 1), 
phantom energy (w < — 1 and quintom (that can across from phantom region to 
quintessence region as evolved) and have time dependent EoS parameter. Some 
other limits obtained from observational results coming from SN la data (Knop 
et al. 35 collaborated with CMBR anisotropy and galaxy clustering statistics 
(Tegmark et al. [35]) are -1.67 < uj < -0.62 and -1.33 < u < -0.79 respec- 
tively. However, it is not at all obligatory to use a constant value of ui. Due 
to lack observational evidence in making a distinction between constant and 
variable u>, usually the equation of state parameter is considered as a constant 
(Kujat et al. [37]; Bartelmann et al. [38] ) with phase wise value — 1,0, — ^ and 
+1 for vacuum fluid, dust fluid, radiation and stiff dominated universe, respec- 
tively. But in general, w is a function of time or redshift (Jimenez [39] ; Das et 
al. [3D]; Ratra and Peebles [H]). For instance, quintessence models involving 
scalar fields give rise to time dependent EoS parameter ui (Turner and White 
[42] ; Caldwell et al. [43]; Liddle and Scherer [44]; Steinhardt et al. [35]). Some 
literature are also available on models with varying fields, such as cosmological 
models with variable equation of state parameter in Kaluza-Klein metric and 
wormholes (Rahaman et al. 46, 47 ). In recent years various form of time de- 
pendent ui have been used for variable A models (Mukhopadhyay et al. [4"5] ; 
Usmani et al. [33]). Recently Ray et al. [49], Akarsu and Kilinc {50] and Yadav 
[511 152] have obtained dark energy models with variable EoS parameter. 

Spatially homogeneous and anisotropic cosmological models play a signifi- 
cant role in the description of large scale behaviour of universe and such models 
have been widely studied in framework of General Relativity in search of a 
realistic picture of the universe in its early stages. Recently, Pradhan et al. 
[53] — [56] and Saha et al. [57] — [59] have studied homogeneous and anisotropic 
B-I space-time in different contexts. In this paper, we have investigated a new 
anisotropic B-I DE model with variable uj by assuming constant deceleration 
parameter. The out line of the paper is as follows: In Section 2, the metric and 
the field equations are described. Section 3 deals with the solutions of the field 
equations in two different cases with physical and geometric behaviour of the 
models. Finally, conclusions are summarized in the last Section 4. 
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2 The Metric and Field Equations 

We consider totally anisotropic Bianchi type-I line element, given by 

ds 2 = -dt 2 + A 2 dx 2 + B 2 dy 2 + C 2 dz 2 , (1) 

where the metric potentials A, B and C are functions of t alone. This ensures 
that the model is spatially homogeneous. 

The simplest generalization of Equation of State (EoS) parameter of perfect 
fluid may be to determine the EoS parameter separately on each spatial axis by 
preserving the diagonal form of the energy momentum tensor in a consistence 
way with the considered metric. Therefore, the energy momentum tensor of 
fluid is taken as 

Tij = diag[T 00 , T n , T 22 , T 33 ] (2) 
Thus, one may parameterize it as follows, 

= diag[p, -p x , -p y , -p z ] = diag[l, -u x , -ui y , -u> z ]p 

= diag[l,-u,-(uj + 6),-(uj + j)]p. (3) 

Here p x ,Py and p z are the pressures, p is the proper energy density and u) x ,u y 
and oj z are the directional EoS parameters along the x, y and z axes respectively. 
ui is the deviation-free EoS parameter of the fluid. We have parameterized the 
deviation from isotropy by setting uj x — uj and then introducing skewness param- 
eters 5 and 7 that are the deviations from u> along the y and z axes respectively. 



The Einstein's field equations (with gravitational units, 8irG = 1 and c = 1) 
read as ^ 

Rij ~ 7^-9 ij — ^^iji (4) 

where the symbols have their usual meaning. In a comoving co-ordinate system, 
Einstein's field equation with Q for B-I metric (JXJ) subsequently lead to 
the following system of equations: 

B C BC 

B + C + BC^-^ (5) 
= -(uj + S)p, (6) 

= -(w + 7)p, (7) 

Here and in what follows an over dot denotes ordinary differentiation with re- 
spect to t. 
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The spatial volume for the model ([T]) is given by 

V 3 = ABC. (9) 

We define a = (ABC) 3 as the average scale factor so that the Hubble's param- 
eter is anisotropic and may be defined as 
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H =a = - 3 [A + B + C)- (10) 
We define the generalized mean Hubble's parameter H as 

H =\{H X +H y + H z ), (11) 

where H x = ^, H y = ^ and H z = are the directional Hubble's parameters 
in the directions of x, y and z respectively. 

An important observational quantity is the deceleration parameter q, which 
is defined as 

ad hoi 

<? = - — ■ ( 12 ) 

The scalar expansion 9, shear scalar a 2 and the average anisotropy parameter 
A m are defined by 

--A 4- 

°°=l(t H t-\ i >°)- <") 

i — 1 

where Ai7, = H t - H(i = 1, 2, 3). 



3 Solutions of the Field Equations 

The field equations ©-(jH]) are a system of four equations with seven unknown 
parameters A, B, C, p, uj, 8 and 7. Three additional constraints relating these 
parameters are required to obtain explicit solutions of the system. 

Firstly we apply the special law of variation for generalized Hubble's parameter 
that yields a constant value of deceleration parameter. Since the line element 
|T]) is completely characterized by Hubble parameter H , therefore, let us con- 
sider that the mean Hubble parameter H is related to the average scale factor 
a by the relation 

H = £a- n =£(ABC)-9, (16) 
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where £(> 0) and n{> 0) are constants. Such type of relations have already 
been considered by Berman [60], Berman and Gomide [61] for solving FRW 
models. Latter on many authors (see, Singh et al. [52] — [M], Singh and Baghel 
[55] . Pradhan and Jotania [55J and references therein) have studied flat FRW 
and Bianchi type models by using the special law for Hubble parameter that 
yields constant value of deceleration parameter. 

From (fTU]) and ([15]). we get 

a = e a - n+1 , (17) 

and 

S= ~e 2 (n~l)a- 2n+1 . (18) 
Substituting HT7J) and (0 into ([11]), we get 

q = n-l. (19) 

We observe that the relation (fl9|) gives g as a constant. The sign of q indicated 
whether the model inflates or not. The positive sign of q i.e. (n > 1) correspond 
to "standard" decelerating model whereas the negative sign of q i.e. < n < 1 
indicates inflation. It is remarkable to mention here that though the current 
observations of SNe la and CMBR favours accelerating models (q < 0), but 
both do not altogether rule out the decelerating ones which are also consistent 
with these observations (see, Vishwakarma [57]). 

Integrating Eq. (fTT|) we obtain the law of average scale factor 'a' as 

a=(n£t + cx)i for n ^ 0, (20) 

and 

a = c 2 e lt for n = 0, (21) 

where c\ and C2 are constants of integration. Thus, the law (|16[) provides 
two types of the expansion in the universe i.e., (i) power-law (|20[) and (ii) 
exponential-law (12T1) . 

Secondly we assume that the component <r\ of the shear tensor a\ is pro- 
portional to the expansion scalar (9), i.e., a\ oc 9. This condition leads to the 
following relation between the metric potentials: 

A=(BC) m , (22) 

where m is a positive constant. The motive behind assuming this condition is 
explained with reference to Thorne |68j , the observations of the velocity-red-shift 
relation for extragalactic sources suggest that Hubble expansion of the universe 
is isotropic today within 30 per cent (Kantowski and Sachs |69j ; Kristian and 
Sachs [75]). To put more precisely, red-shift studies place the limit 

^<0.3, 
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on the ratio of shear a to Hubble constant H in the neighbourhood of our 
Galaxy today. Collins et al. [71] have pointed out that for spatially homoge- 
neous metric, the normal congruence to the homogeneous expansion satisfies 
that the condition % is constant. 

Thirdly we assume that the deviations from u> along y and z axes are same 
i.e. 7 = 5. 

Subtracting (j6|) from ((7]), and taking integral of the resulting equation two times, 
we get 

L c A J (ABC)~ 1 dt , (23) 



— = c 3 exp 



where C3 and c 4 are constants of integration. 

In the following Subsections (3.1) & (3.2), we consider two types of solutions, 
i.e., for power-law and exponential form volumetric expansion. 



3.1 Case (i): When n^O i.e. model for power-law expan- 
sion 

After solving the field equations (J6j) - ([8]) for the power-law volumetric expansion 
([20]) by considering Eqs. (|22|) and (|23|) . we obtain the scale factors as follows 



A(t) = (n£t + c 1 )^ m + 1 i , 



B (t) = y/c3(nit + ci) 2 "C»+« exp 



1 a 
C(t) = — — (nft + ciW m + 1 ) exp 
VC3 



2i{n - 3) 
c 4 



2£(n - 3) 
provided n 7^ 3. 

Hence the model (JTJ) is reduced to 

ds 2 = -df+init+a) dx 2 +c 3 (nlt+ci) exp 



{nlt + d)- 
(nit + c 1 ) 1 



(24) 
(25) 
(26) 



c 4 



£(n - 3) 



dy 2 



1 3 
-I (nit + ci)"( m + 1 > exp 

C3 



c 4 



-(n& + ci) j 



£(n - 3) 

The rate of expansion if^ in the direction of x, y and z read as 
fl* = - = rCnft + ci) -1 , 



dz 2 



j4 TO + 1 



3/' 



JJ 

v ~ B ~ 2(to+ 1) 



1 



(rdt + ci) 1 + -c 4 (n£t + c x ) ~, 



(27) 

(28) 
(29) 
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C 3£ 1 
U z(m + lj 2 



(30) 



Thus the Hubble's parameter H, scalar of expansion 9, shear scalar a and the 
average anisotropy parameter A m are given by 



a 2 = 



6 = 3H = U{nlt + ci)- x , 

U 2 (2m- 1) 2 , . . 2 1 o, „ 
2(m + l) 2 {ntt + Cl) + 4°^ + ClY 



A — 



(2m- 1) 



(m+1) 2 6^ 2 



+ -^(nft + Cl ) ; 



(31) 
(32) 

(33) 
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Figure 1: The plot of EoS parameter u> versus t in power- law expansion for 
q < 0. Here t = 0.33, n = 0.5, m = 2, c 4 = 2 



The energy density of the fluid can be find by using Eqs. (|24l) - (1261) in 

9£ 2 (4m+ 1) 



/'= 4 ( TO+ ^2 ' (n£t + a) 2 - -cl(n£t + Cl y 



(34) 
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Figure 2: The plot of EoS parameter ui versus t in power-law expansion for 
q > 0. Here £ = 0.33, n = 2, m = 2, c 4 = 2 



Using Eqs. (j25|) . f[26]) and ([34| in ([5j> , the EoS parameter ui is obtained as 



U 2 



-4n(ro+l)+9 



4(m+l) 



(n& + ci)- 2 + ±cf (n# + ci)-« 



\c\{nit + ex)"* - ^^^(nft + cx)-a 



(35) 



Using Eqs. (|24 |) - (|26|) . (|34|) and (|35|) in either (|6J) or (J7J)j the skewness parameters 
S and 7 (i.e. deviation from uj along y-axis and z-axis) are computed as 



1„2 
4 



(36) 



From equation (|35[) . it is observed that the equation of state parameter ui is 
time dependent, it can be function of redshift z or scale factor V as well. The 
redshift dependence of uj can be linear like 



u)(z) = uio + uj 1 z, 



(37) 
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Figure 3: The plot of energy density p versus t in power-law expansion for q < 0. 
Here £ = 0.33, n = 0.5, m = 2, c 4 = 2 



Weller and Albrecht 



(38) 



with w 1 = (^)z = (see Refs. Huterer and Turner 
[30 ) or nonlinear as 

w{z) = UJ + — — , 
1 + z 

(Polarski and Chavellier [72]; Linder [73]). The SNe la data suggests that 

— 1.67 < u) < —0.62 (Knop et al. [3S]) while the limit imposed on a; by a combi- 
nation of SNe la data (with CMB anisotropy) and galaxy clustering statistics is 

— 1.33 < cj < —0.79 (Tegmark et al. [36 ). So, if the present work is compared 
with experimental results mentioned above, then one can conclude that the limit 
of u) provided by equation (|35[) may accommodated with the acceptable range 
of EoS parameter. Also it is observed that at t = t c , uj vanishes, where t c is a 
critical time given by 



tr 



1 f £y/3[4n(m + l)-9]\ n c i 
ni \ cAm + 1) / nf 



Thus for this particular time our model represents a dusty universe. 



(39) 
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Figure 4: The plot of energy density p versus t in power-law expansion for q > 0. 
Here £ = 0.33, n = 2, m = 2, c 4 = 2 



For the value of u> to be in consistent with observation (Knop et al. [35]). 
we have the following general condition 



where 



and 



h < t < t 2 , 



h --£ 



h -~£ 



1.121 



y/4n(m + 1) + 5(4m + 1) 
c 4 (to + 1) 



y/4n(m + 1) + 7.44m - 7.14 
c 4 (m + 1) 



£i 



ni' 



(40) 
(41) 

(42) 



For this constrain we obtain —1.67 < uj < —0.62 which is in good agreement 
with the limit obtained from observational results coming from SNe la data 
(Knop et al. [35] ). We note that the earlier real matter at t < t c , where uj > 
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Figure 5: The plot of cosmological constant A versus t in power-law expansion 
for q < 0. Here £ = 0.33, n = 0.5, m = 2, c 4 = 2 



later on at t > t c , where lu < converted to the dark energy dominated phase 
of universe. 

We also observe that if 



to 



1 



1.221 



^An(m + 1) + Q(2m- 1) 



04(771 + 1) 



72i' 



(43) 



for i — to, u = — 1 (i.e. cosmological constant dominated universe), and when 
t < to, lu > — 1 (i.e. quintessence), and for t > to, to < —1 (i.e. super 
quintessence or phantom fluid dominated universe) (Caldwell |74j ) . 



Figures 1 & 2 depict the variation of equation of state parameter (lu) ver- 
sus cosmic time (t) in the two modes (i.e. q < and q > 0) of evolution of 
the universe, as a representative case with appropriate choice of constants of 
integration and other physical parameters using reasonably well known situa- 
tions. From Figure 1 & 2, we conclude that in early stage of evolution of the 



12 



120 



100 



80 



60' 



40' 



20 




Figure 6: The plot of cosmological constant A versus t in power-law expansion 
for q > 0. Here I = 0.33, n = 2, m = 2, c 4 = 2 



universe, the EoS parameter was positive (i.e. the universe was matter dom- 
inated) and at late time it is evolving with negative value (i.e. at the present 
time). The earlier real matter later on converted to the dark energy dominated 
phase of the universe in both accelerating & decelerating modes. 



From Eq. (|34l) . we note that energy density of the fluid p{t) is a decreasing 
function of time and p > when 



t > 



( C4 (m 



1 



V 3£y/4m + 1 



ci 



(44) 



Figures 3 &: 4 are the plots of energy density of the fluid (p) versus time in 
accelerating and decelerating modes of the universe respectively. In both modes 
we observe that in early universe p is negative and suddenly reaches to maxi- 
mum positive value and then decreases with time and at late time it approaches 
to zero. 
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In absence of any curvature, matter energy density fi m and dark energy Q\ 
are related by the equation 

Q m + Q A = 1, (45) 
where fl m = and fl\ = Jk?. Thus equation (1431) . reduces to 



3H* " A 3fP> 

P A 



3H 2 3H 2 



1. (46) 



Using equations ([31]) and in equation (|4"6")l . the cosmological constant is 

computed and obtained as 

3£ 2 (4m 2 - 4m+ 1) , „ . 2 1 „ , 6 
A = 4(m+l) 2 ( ^ + Cl) + 4^ + Cl) " • (47) 

From eq. (I47|) we observe that A is a decreasing function of time and it is always 
positive for m > |. 

In recent time the A-term has interested theoreticians and observers for var- 
ious reasons. The nontrivial role of the vacuum in the early universe generate 
a A-term that leads to inflationary phase. Observationally, this term provides 
an additional parameter to accommodate conflicting data on the values of the 
Hubble constant, the deceleration parameter, the density parameter and the age 
of the universe (for example, see the references Gunn and Tinsley [75]; Wampler 
and Burke 76 ) . The behaviour of the universe in this model will be determined 
by the cosmological term A, this term has the same effect as a uniform mass 
density p e // = — A which is constant in time. A positive value of A corresponds 
to a negative effective mass density (repulsion). Hence, we expect that in the 
universe with a positive value of A the expansion will tend to accelerate whereas 
in the universe with negative value of A the expansion will slow down, stop and 
reverse. In a universe with both matter and vacuum energy, there is a com- 
petition between the tendency of A to cause acceleration and the tendency of 
matter to cause deceleration with the ultimate fate of the universe depending 
on the precise amounts of each component. This continues to be true in the 
presence of spatial curvature, and with a nonzero cosmological constant it is no 
longer true that the negatively curved ("open") universes expand indefinitely 
while positively curved ( "closed" ) universes will necessarily re-collapse - each of 
the four combinations of negative or positive curvature and eternal expansion or 
eventual re-collapse become possible for appropriate values of the parameters. 
There may even be a delicate balance, in which the competition between matter 
and vacuum energy is needed drawn and the universe is static (non expanding) . 
The search for such a solution was Einstein's original motivation for introducing 
the cosmological constant. 

Figures 5 & 6 are the plots of cosmological constant A versus time in ac- 
celerating and decelerating mode of the universe respectively. In both mode 
we observe that cosmological parameter is decreasing function of time and it 
approaches a small positive value at late time (i.e. at present epoch). Recent 
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cosmological observations (Garnavich et al. Q][2]; Perlmutter et al. [3]— [5]; 
Riess et al. [6] [77l [78]; Schmidt et al. [7]) suggest the existence of a positive 
cosmological constant A with the magnitude A(GTi/c 3 ) « 10 -123 . These ob- 
servations on magnitude and red-shift of type la supernova suggest that our 
universe may be an accelerating one with induced cosmological density through 
the cosmological A-term. But this does not rule out the decelerating ones which 
are also consistent with these observations (Vishwakarma [67]). Thus the nature 
of A in our derived DE model is supported by recent observations. 

From Eqs. (|28p - (|3"3")) . it can be seen that the spatial volume is zero at t = — 
and it increases with the cosmic time. The parameters Hi, H, 9 and a diverge 
at the initial singularity. There is a Point Type singularity (MacCallum [79j ) 
at t — — ^ in the model. The mean anisotropic parameter is an increasing 
function of time for n > 3 whereas for n < 3 it decreases with time. Thus, the 
dynamics of the mean anisotropy parameter depends on the value of n. Since 
§2- = constant (from early to late time) , the model does not approach isotropy 
through the whole evolution of the universe. 



3.2 Case (ii): When n = i.e model for exponential ex- 
pansion 

After solving the field equations ([6]) - ((HJ for the power-law volumetric expansion 
(p?Tj) by considering Eqs. ([22)) and (|23l) . we obtain the scale factors as follows 



A(t) = c 2 m+1 exp 



3m£ 



B(t) - ^2 exp 



to + 1 

31 



t-^ e - ut 



C(t) 



■ exp 



2(m + 1) Uc\ 



2 (to + 1) Uc\ 



Hence the model {U is reduced to 

6m£ 



-dt +C2 exp 



to + 1 



t I dx +C3C2 exp 



3^ ^ _ e 4 -set 



(m + l)~ Uc\ 



(48) 
(49) 
(50) 

dy 2 



C3 



• exp 



(to + 1) Uc\ 



dz l 



(51) 



The rate of expansion Hi in the direction of x, y and z and the Hubble 
parameter are obtained as 

Hx = — , (52) 



TO + 1 



H,i 



3t 



C 4 -3lt 



2(to+1) 2c?, 



(53) 
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Figure 7: The plot of EoS parameter uj versus t in exponentialr-law expansion 
for q = -10. Here t = 0.33, n = 0, m = 2, c 4 = 2 



£T, = ? ^ . - ^e~ m . (54) 
2(m + 1) 2c\ 



Thus the Hubble's parameter H, scalar of expansion 6, shear scalar a and the 
average anisotropy parameter A m are given by 

H = £, (55) 

6 = U, (56) 

(57) 



3^ 2 (2m-l) 2 4 



o 2 = — 



4(m+l) 2 4e| 



(2m- l) 2 



2 
2 

4 „-e 



m 2(m+l) 2 6^ 2 c| 
Using Eqs. (|4"5)) - (f5T))) in JSJ, the energy density of the fluid is obtained as 



(58) 



9^ 2 (4m + 1) c 2 
4(m + l) 2 4cSj 



e" M *. (59) 
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Using Eqs. P^|) . ([50]) and (|5§f in (J3J) , the EoS parameter cj is obtained as 

27£ 2 i e 2 6ft 
4(m+l)" + i^f e 

w = —5 r ■ (60) 

c l „-6ft _ 9l 2 (4m+l) V ; 

iSJ e 4(m+l)^ 

Using (|i8 |) -(|50" ]l and ([59]). ([60]) in either § or 0, the skewness parameters 6 
and 7 are computed as 

9l 2 (2TO 2 +m-l) 

S = 7 = 2(m+1) ' (61) 

7 9£ 2 (4m+l) _ c 2 6ft ■ ^ 
4(m+l) 2 23> c 

It is observed that for t = i c , a; vanishes, where i c is a critical time given by 

1 / c!(m + l) 2 \ 
* C ~6^H 27c^ J' ^ 
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10 



Figure 9: The plot of cosmological constant A versus t in exponential-law ex- 
pansion for q = — 1. Here I = 0.33, n — 0, m = 2, C4 = 2 



Thus for this particular time our model represent a dusty universe. 



For the value of u to be in consistent with observation (Knop et al. 
we have the following general condition 



where 



and 



ti < t < t 2 , 



t 4(m + lf 

1 U V3-37c^ 2 (6.7m- 1.33) 



t 2 = — In 
U 



cl{m + lf 



5.6c^ 2 (2.4m- 2.38) 



[35]), 



(63) 



(64) 



(65) 



For this constrain we obtain —1.67 < u> < —0.62, which is in good agreement 
with the limit obtained from observational results coming from SNe la data 
(Knop et al. [35]). We note that the earlier real matter at t < t C) where oj > 0, 
later on at t > t c , where u) < 0, converted to the dark energy dominated phase 
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of universe. 



It is also observed that if 

1 ( 4(m + l) 2 \ 
^^{ toZPim-l) )' (66) 

for t = t$, lu = —I (i.e. cosmological constant dominated universe), and when 
t < to, lj > — 1 (i.e. quintessence), and for t > t , uj < — 1 (i.e. super 
quintessence or phantom fluid dominated universe) (Caldwell [74]). 

The variation of equation of state parameter (w) with cosmic time (t) is clearly 
shown in Figures 7, as a representative case with appropriate choice of con- 
stants of integration and other physical parameters using reasonably well known 
situations. From Figure 7, we conclude that in early stage of evolution of the 
universe, the EoS parameter u was positive (i.e. the universe was matter dom- 
inated) and at late time it is evolving with negative value (i.e. at the present 
time). The earlier real matter later on converted to the dark energy dominated 
phase of the universe. 

From (15^1) . it is observed that p(t) is a slowly increasing function of time for 
large t. When t — > oo, p becomes a constant. The energy condition p > is 
satisfied under the condition 

> C4(TO + 1) (67) 
3£4^(4m+l) 

This behavior of p(t) is clearly depicted in Figure 8. This means there is no 
density evolution in this set up. A possible reason for no evolution of energy 
density is that the universe could be much rapid in which the matter do not get 
time to readjust on expanding range or other unknown dominated effect which 
not incorporated in potential functions of this space-time for exponential-law. 

Using equations (|55|) and (j59")l . in equation PB^I . the cosmological constant is 
obtained as 

3l 2 (4m 2 -4m + l) c\ m 
A - 4(m+l)2 + 44 6 ■ (68) 

From Eq. we observe that A is a decreasing function of time and it is 

always positive for m > i . This behavior of cosmological constant A is shown 
in Figure 9. Thus, in this case also we observe that A is consistent with recent 
cosmological observations (Garnavich et al. [HH]; Perlmutter et al. [3]— [5]; 
Riess et al. [51173 ES]; Schmidt et al. 0) 

From above results, it is observed that the physical and kinematic quantities 
are all constant at t = 0. From Eqs. (|52p - ([55]) . we observe that the kinematic 
parameters tend to zero as t —¥ oo. The expansion in the model is uniform 
throughout the time of evolution. Since ^ = constant (from early to late 
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time), the model does not approach isotropy at any time. The derived model is 
non-singular. 

4 Concluding Remarks 

Two new anisotropic B-I DE models with variable EoS parameter uj have been 
investigated which are different from the other author's solutions. In the derived 
models, uj is obtained as time varying which is consistent with recent observa- 
tions (Knop et al. [35]; Tegmark et al. [35]). The proposed law of variation for 
the Hubble's parameter yields a constant value of deceleration parameter. The 
law of variation for Hubble's parameter defined by (fT()]) for B-I space-time gives 
two types of cosmologies, (i) first form for (n ^ 0) gives the solution for positive 
value of deceleration parameter for n > 1 and also negative value of decelera- 
tion parameter for < n < 1 indicating the power law expansion of the universe 
whereas (ii) second one (for n = 0) gives the solution for negative value of de- 
celeration parameter, which shows the exponential expansion of the universe. 
The power law solution represents the singular model where the spatial scale 
factors and volume scalar vanish at t = — All the physical parameters are 
infinite at this initial epoch and tend to zero as t — > oo. There is a Point Type 
singularity (MacCallum [79] ) at t = — in the first model. The exponential 
solutions represents singularity free model of the universe. In both cases, since 
= constant (from early to late time), the models do not approach isotropy 
through the whole evolution of the universe. 

The main features of the models are as follows: 

• The DE models are based on exact solutions of the Einstein's field equa- 
tions for the anisotropic B-I space-time filled with perfect fluid with vari- 
able EoS parameter uj. The exact solutions of the Einstein field equations 
have been obtained by assuming two different volumetric expansion laws in 
a way to cover all possible expansions: namely power-law and exponential- 
law expansion. To my knowledge, the literature has hardly witnessed this 
sort of exact solutions for the anisotropic B-I space-time. So the derived 
DE models add one more feather to the literature. 

• In both cases (i) & (ii), it is observed that, in early stage, the EoS pa- 
rameter uj is positive i.e. the universe was matter dominated in early 
stage but in late time, the universe is evolving with negative values i.e. 
the present epoch (see, Figures 1 & 2). Thus our DE models represent 
realistic models. 

• In both cases the DE models present the dynamics of EoS parameter ui 
provided by Eqs. (|35p and (|60l) respectively may accommodated with the 
acceptable range —1.67 < ui < —0.62 of SNe la data (Knop et al. [35]). It 
is already observed and shown in previous sections that for different cosmic 
times, we obtain cosmological constant dominated universe, quintessence 
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and phantom fluid dominated universe (Caldwell [73]). representing the 
different phases of the universe through out the evolving process. Unlike 
Robertson- Walker (RW) metric Bianchi type metrics can admit a DE 
that wields an anisotropic EoS parameter according to the characteristics. 
The cosmological data - from the large-scale structures [SU] and Type 
la supernovae [75J HI] observations - do not rule out the possibility of 
an anisotropic DE either [82] 155] . Therefore, one can not rule out the 
possibility of anisotropic nature of DE at least in the frame-work of B-I 
space-time. 

• Our DE models are of interest because in both cases the nature of de- 
caying vacuum energy density A(t) are supported by recent cosmological 
observations (Garnavich et al. [T] Perlmutter et al. [3J— [5]; Riess et 
al. jBH77J[75J; Schmidt et al. [7 ). These observations on magnitude and 
red-shift of type la supernova suggest that our universe may be an accel- 
erating one with induced cosmological density through the cosmological 
A-term. 

• Though there are many suspects (candidates) such as cosmological con- 
stant, vacuum energy, scalar field, brane world, cosmological nuclear- 
energy, etc. as reported in the vast literature for DE, the proposed model 
in this paper at least presents a new candidate (EoS parameter) as a 
possible suspect for the DE. 
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